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Abstract. The anomalies of a very general class of non local Dirac operators are com- 
puted using the ('-function definition of the fermionic determinant and an asymmetric 
version of the Wigner transformation. For the axial anomaly all new terms introduced 
0^ ' by the non locality can be brought to the standard minimal Bardeen's form. Some 

extensions of the present techniques are also commented. 
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^ INTRODUCTION 

(N' 

2 ! Local field theories provide the commonly accepted setup where the implementa- 

Q\ [ tion of space-time symmetries becomes rather simple. On the other hand, effective 

^ ■ theories are not necessarily local, although an appropriate choice of degrees of free- 

dom can make them almost local [1]. In low energy QCD, light quarks and gluons 
, are dressed by the interaction in a way that the effective theory looks highly nonlo- 

Q I cal [2], and a kind of dynamical perturbation theory would be needed [3]. In terms 

1^ • of pions and (heavy) nucleons the theory becomes weakly nonlocal and a chiral 

perturbation theory becomes of practical interest [4]. In a Dyson-Schwinger setting 
k>( I [5] most information about such a non local theory comes from the constraints im- 

^ I posed by the relevant Ward and Slavnov- Taylor identities [5], perturbation theory 

to some finite order and hadronic phenomenology. These approaches are necessary 
if one wants to know, for instance, about the momentum distribution of a quark in 
a hadron; neither perturbative QCD nor chiral perturbation theory can properly 
handle this problem. In this regard, anomalies are particulary interesting because 
their existence is linked to a violation of classical symmetries by high energy regu- 
lators although their physical effect is formulated as a low energy theorem. At the 
one loop level, anomalies for nonlocal models have been previously discussed [6-8] 
for some specific processes like e.g. vr" -^ 27, 7 — i> Svr and 2K —* Svr. We refer 
to those works and [9] for further motivation. Rather than computing all specific 
processes one by one we prove that the new terms generated by the non locality 
can be subtracted by adding suitable counterterms. 
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I THE ONE LOOP EFFECTIVE ACTION 

Our starting point is the effective action of Dirac fermions in tlie fiat Euclidean 
space-time endowed with internal degrees of freedom collectively referred to as 
"flavor", 

W{D) = -log fv^Vijexpl- f d^x4){x)D4j{x)\ = -TrhgD. (1) 

Here, Tr stands for trace over all degrees of freedom and D is the Dirac operator to 
be specified below. The definition of the fermion determinant requires some renor- 
malization of the ultraviolet divergences. The (consistent) chiral anomaly is defined 
as the variation of the effective action under infinitesimal chiral transformations, 
given by 

tpix) -^ e^^-^"^5^(x) , tpix) -^ i,{x)e-'^-'"^<'^ (2) 

where a{x) and j3{x) are Hermitian matrices in fiavor space only, regarded as multi- 
plicative operators on the fermionic wave functions. The particular cases a = and 
/9 = correspond to vector and axial transformations, respectively. This induces 
transformations of the Dirac operator 

which infinitesimally become 

6D = 6vD + 6aD = [i(3, D] - {laj^, D} . (4) 

Since we will be considering a C- function renormalization of W (see below), there 
will be no vector anomaly, 

5yW = 0, 6aW = Aa . (5) 

Correspondingly, the same current conservation formulas valid for the local case 
can be written here, 



= J d''x{^{x)[tf3,D]^{x))Q (6) 

-Aa = f d^x {^{x){ia-f5, D}^{x))q . (7) 



(The symbol ( )q stands for quantum vacuum expectation value.) In particular 
the term 7^-P/i in D in the right-hand side yields, after integration by parts, the 
divergence of the fermionic vector and axial currents whereas the other terms in 
D, local and non local, represent the explicit chiral symmetry breaking due to the 
external fields. On the other hand, the left-hand side shows the anomalous breaking 
of the axial current conservation. 



The class of Dirac operators to be considered here is 

D = Dl + M. (8) 

The term Di, the local component of D, is a standard Dirac operator 

Dl = 7m^m + Y (9) 

We will follow the conventions of [10,11]), P^ = id^ and Y is an arbitrary matrix- 
valued function in flavor and Dirac spaces. Y is a function of the position operators 
X^, deflned by Xfj_ip{x) = x^ip{x), so that "K is a multiplicative operator in the 
Hilbert space of fermions. The term M is a purely non local ^ , more precisely 
bilocal, operator also with arbitrary structure in flavor and Dirac spaces, 

{Ytl)){x) = Y{x)iIj{x) , {MiIj){x) = f d^yM{x,y)ilj{y) . (10) 



More restrictive assumptions on M are spelled out in ref. [11]. For the purpose of 
doing detailed calculations we will assume that the non local operator M admits 
an expansion in inverse powers of P^ for large P^ of the form 

P P P P P P 

M = M^^ + M^,^ + M^,, ^f' + • • • (11) 

The coefficients Ai"^j...^^ are multiplicative operators and they are completely sym- 
metric under permutation of indices. For convenience the P^ has been put at the 
right ^ . The better way to obtain the transformation properties of Nl is by intro- 
ducing a family of operators associated to M as 

M[p) = e'P^Me~'P^ (12) 

where the momentum p^ is just a constant c-number. Effectively, M{p) corresponds 
to make the replacement Pfj, ^ Pfj. + p^ in M. The function M{p) admits an 
expansion in inverse powers of p^ similar to that in eq. (11), namely 

M{p) = M,^ + M,,^ + ■■■. (13) 

We will adopt the ^-function renormalization prescription combined with an 
asymmetric Wigner transformation. This method, as well as several of its ap- 
plications, is presented in great detail in [10]. The ^-function effective action is 
given by [13,14] 



'^' M is softer in the ultraviolet sector than any multiplicative operator, that is, the distribution 
M{x,y) is less singular than the Dirac delta 5{x — y). 

■^^ This choice does not exhaust all possible non local operators, but it is realistic since it ac- 
commodates the operator product expansion estimate of the quark self-energy S(p^) ^p^^oo 
{\ogp'^Y~^ / p^ with d the anomalous dimension of the quark condensate V'V' (see ref. [12]) 



© 


X 




— 


r V_ 

X 


^ X 



FIGURE 1. The F contour in the complex z-plane. Crosses represent isolated eigenvalues of D. 



WiD) = -Trlogl? = -—Tt(D' 

as 



's=0 ' 



(14) 



where s = is to be understood as an analytical extension on s from the ultraviolet 
convergent region Re(s) < —D. The operator D^ can be obtained from 

dz z^ 



D' = - 



r 27ii D — z 



(15) 



where the integration path F starts at — oo, follows the real negative axis, encircles 
the origin z = clockwise and goes back to — oo (see figure 1). The key point 
is that for sufficiently negative s there are no ultraviolet divergences and formal 
operations become justified. By construction, the ^-function renormalized effec- 
tive action is invariant under all symmetry transformations associated to similarity 
transformations of D, thus in particular it is vector gauge invariant. On the other 
hand, the variable z plays the role of a mass and hence breaks explicitly both chiral 
and scale invariance. 



The operator {D- 



\-i 



can be conveniently expressed by means of an asymmetric 



version of the Wigner representation [15]. For any operator A, let 

{x\A\p) 



A{x,p) = f d^ye'yP{x\A\x-y) 



{x\p) 



(16) 



be its (asymmetric) Wigner representation. \p) is the momentum eigenstate with 
{x\p) = e~*^^. From this definition 



x\A\x) = j -——^A{x,p), TrA = y 



d^xd^p 



D 



ti A{x,p) 



(17) 



where tr acts on internal and Dirac spinor degrees of freedom only, and the product 
of two operators satisfy the following formula 



{AB){x,p) = expiid^ ■ d^)A{x,p)B{x,p) 



'It 



where dp acts only on the p-dependence in A(x,p) and d^ on the x-dependence in 
B{x,p). Let the propagator or resolvent of D, be G{z) = {D — z)~^ and G{x,p; z) 
its Wigner representation. Applying Eq.l8 to {D — z)G = 1 one obtains 



G{x,p-z) = {x\{i)+D-z)-^\{)) (19) 

where |0) is the state of zero momentum, (a;|0) = 1. In practice this imphes that 
idfj, derivates every x dependence at is right, until it annihilates |0), iS^lO) = 0. 
This method is very efficient for a derivative expansion for it computes directly, 
that is non recursively, each of the terms [10]. 

The definition given for A{x,p) is not gauge covariant because \p) is not. We will 
consider only local objects of the form {x\f{D)\x) as given by the formula 

{x\f{D)\-) = -j£^i ^/WG(-'^; -) (3-7) 

We will assume that the function f{z) is sufficiently convergent at infinity or else 
that it can be obtained as a suitable analytical extrapolation from a parametric 
family f{z, s) in the variable s. In either case the integration over z should be 
performed in the first place, to yield the Wigner representation of the operator 
f{D). Afterwards, the p integration is carried out, corresponding to take the 
diagonal matrix elements of {y\f{D)\x), hence restoring gauge covariance. This 
obviously means that the gauge breaking piece of G{x,p; z) is a total derivative in 
the momentum variable. Recently, a method has been developped [16] where this 
total derivative is, by construction, gauge invariant. 



II ANOMALIES 

Because the chiral transformations are local, both D^ and M transform covari- 
antly separately, that is, 

6Dl = [tP, Dl] - {^a75, Dl} , 6M = [z/3, M] - {^75, M} . (20) 

Note that the bilocal structure of M implies that local factors at each 
side of the operator are taken at different points, i.e. M{x, x') -^ 

gi/3(a;)-ia(a;)75 ]\^(rj. ^'\^-il3{x')-ia{x')^5 

The two lowest coefficients are given by 

M^ = M^, M^, = M^, + t^^p^MpP, , (21) 

where we have introduced t^yptj = S^^Sp^j — S^pSy^j — S^^jS^p. It should be noted 
that the coefficients Af^^...^^ are not multiplicative operators. One useful property 
of M{p) is that it transforms covariantly under chiral transformations. Indeed, if 
Mfi = Q1MQ2 for two multiplicative operators ^1,2, 

Mn{p) = e*^^Mf,e-*P^ = Q,M{p)Q2 ■ (22) 

As a consequence, the coefficients are also chiral covariant 



SM,,...,„ = [iP, M,,...,„] - {^aJ5, M,„„,J . (23) 

From here it is immediate to derive the transformation of the original coefficients 
-^/ii.../i„- For the two lowest order coefficients one finds 

6M^ = [i(3, M^] - {2a75, M^] , (24) 

5M^, = [if3, M^,] - {^75, M^,} + t^,^„Mp{d^(3 + d^a-/^) . 

In general, the variation of each coefficient involves those of lower order. 

The scale transformation tp{x) -^ e^"'^^^^^^^'^^{e~°'^ x) induces the corresponding 
transformation in D, namely, 

Y{x) ^ e-°^-r(e-"^a:) , M,„„,^{x) ^ e-"^("+^)M^,...^„(e-"^x) . (25) 

Infinitesimally, it implies 

6sD = -as{D - t[X,P„ D]) . (26) 



A Chiral anomaly 

Due to the ( regularization the chiral anomaly becomes an axial anomaly, 

AA = TT{2ta^,D-')^^^. (27) 

Using the Wigner transformation technique [10], the anomaly can be written as (a 
similar expression holds for the effective action) 

Aa = - /t^ / —z'tT(0\2ta^5^^ 10) . (28) 

Here tr stands for trace over Dirac and flavor degrees of freedom, |0) is the zero 
momentum state normalized as (a:|0) = 1, thus -P^|0) = (0|P^ = 0. Further 

D(p) = e'P^De"P^ =i> +Dl + M{p) . (29) 

The integration over z should be performed first, since it defines the operator D^, 
then the integral over p which corresponds to take the trace over space-time degrees 
of freedom and finally s is to be analytically extended to s = 0. The simplest way 
to proceed is to introduce a mass term, i.e., to apply the formula to the Dirac 
operator D + m and then make an expansion in powers of D^ + M{p), letting 
m — > at the end. In this way the following expression is derived 

^ f d^p f dz ^ ,^,^. {DL + M{p))U^+z-m){DL + M{p))) 
^^ = ]^_J (2^1 2^^ *^^°|2^"^^ ipKiz-my)^-^ -^'^ 

(30) 



s=0,m=0 



Simplification has been achieved by using the cychc property for the trace in Dirac 
space. 

After an angular average over p^, the indicated integrals on p^ and z can be 
carried out directly with the integral Ji given in [10]. The result for the four 
dimensional chiral anomaly is 

Aa = (2^75 [^Dl + -Dl{^^, DlYDl + -Ml + -{D^, M^^ 

+ i (m^{7m, Dl]Dl + £>l{7m, M^]Dl + £>l{7m> Dl]M,) ] ) (31) 

The notation (/) stands for (/) = /^^|g/2 tr(0|/(X)|0) . Note that, even for non lo- 
cal Dirac operators, the anomaly is a local polynomial of dimension D constructed 
with P^ and the external fields Y and TVf^j...^^. This is a general property of all 
anomalies since only ultraviolet divergent terms can contribute to them. The ex- 
pressions found for the anomaly can be put in a more usual form, in terms of vector 
and axial fields, scalar fields, etc "^, but it is preferable to use compact notation. 
Since the regularization preserves vector gauge invariance, the axial anomaly is also 
invariant. In our expression for the anomaly, this is a direct consequence of the 
operators there being multiplicative. Indeed, any operator / constructed with the 
gauge covariant blocks Dl and Af^^...^^ is also covariant, i.e., / -^ VlfVt^^ . If in 
addition / is multiplicative (/) is invariant. Note that ( ) is not a trace and so the 
cyclic property does not hold for arbitrary non multiplicative operators. 

B Trace anomaly 

The corresponding trace anomaly, within the ^-function method is [10] 

As = 6sW = asTT{D')s=^. (32) 

The calculation is entirely similar to that of the axial anomaly, yielding 



2 ^ 12 



As = as{ ^D\ + -(£>i{7M, DlY + {7^, Dl}DI{^„ D^} + {7^, DlYdD 



+ 1({7^, £)l}^{7., DlY + ({7^, £)i}{7., Dl})' + {1,, £>l}{7., i^L}'{7^, Dl}) 

+ ^(7^M^Di + ^^DlM^ + -i^DlM^Dl) 

- — (M^7^Di + Dl^^M^ + Dl{i^, M^}Dl + DlI^DlM^ + M^DlI^^Dl) 

+S,iual3{—{l,iM^-fo^DL'yi3DL + ^1 ^D l1 aM ^^i pD l + -f^DL-faDLlf^Mu) 
1 _ _ 1 _ 1 _ 

+— 7^M,7„M^ + — {7„M^„ ^pDl} + —^^M.afs) ) • (33) 



'^' See ref [10] for an explicit expression in the local case and the remarks of ref [11] for the non 
local case. 



Where S^^^a/s = ^fiu^a/s + S^laSu/3 + ^^ip^av The result is again a local polynomial of 
dimension D in the external fields and their derivatives. Unlike the axial case, the 
coefficients Mfj,ua in four dimensions do contribute to the scale anomaly. 

Because scale and chiral transformations commute (in a properly defined sense), 
the crossed variations SsAv,a and 6v,A'A.s coincide and they vanish since the axial 
anomaly is scale invariant. Thus the scale anomaly must be chiral invariant. The 
vector gauge invariance of the previous expressions is easy to check noting that the 
operators inside ( ) are multiplicative. Axial invariance is much more involved in 
general. In four dimensions it is relatively easy to check that the trace anomaly is 
axially invariant in the particular case of M^ = 0, which defines a class of operators 
invariant under chiral and scale transformations. 



C Counterterms and minimal form of the anomaly 

Presumably due to its topological connection [17], the axial anomaly is a very 
robust quantity. It is not affected by higher order radiative corrections [18], and 
remains unchanged at ffnite temperature and density [19]. It gets no contributions 
from scalar and pseudo scalar ffelds [20], tensor ffelds [21,22,10] or internal gauge 
ffelds, i.e, transforming homogeneously under gauge transformations [23,24]. In 
all known cases, the anomaly only affects the imaginary part of the effective action 
in Euclidean space and only involves vector and axial fields. The counter terms 
can always be chosen so that the axial anomaly adopts the minimal or Bardeen's 
form [20]. Not surprisingly, the new terms introduced in the anomaly by the non 
local component of the Dirac operator are also unessential, that is, they can be 
removed by adding a suitable local and polynomial counter term to the effective 
action. In other words, all new terms can be derived as the axial variation of an 
action which is a polynomial constructed with the external fields Y and A^f^^ ^^ 
and their derivatives. The dimension of the polynomial can be at most D. 

The general proof that the anomaly can always be brought to Bardeen's form 
has already been presented [11]. The actual construction of the counter terms can 
be done using the method in ref. [10] (see some further details in ref. [27]). One 
interesting insight in the local case [10] is that the needed counterterms require 
not only the Dirac operator D but also its adjoint D^ not related to the original 
theory. 

The scale anomaly is already minimal. It can be modified by adding polynomial 
counter terms of dimension smaller than D but this would add terms of the same 
type to the scale anomaly. 

Ill EXTENSION TO FINITE TEMPERATURE 

The Wigner transformation method combined with the ^-function regularization 
has been further extended to the finite temperature case in [28,29]. As is well 
known, in the imaginary time formulation of finite temperature field theory, the 



field configurations are periodic or antiperiodic functions of tlie Euclidean time for 
bosons and fermions respectively and thus the frequency running in the fermion loop 
takes discrete values only, cj„ = 7r(2n+l)T (where T stands for the temperature and 
n is any integer) which are known as Matsubara frequencies. At finite temperature, 
the trace of an operator f{xo,x;ido,i'V), acting on the Hilbert space of rf + 1- 
dimensional fermions with possible internal degrees of freedom, becomes 

Trif) =TJ2 f ^^tT{0\f{xo,x;u^ + tdo,p + tV)\0) . (34) 

This formula generalizes that for zero temperature. Note that / is a periodic 
function of xq, and |0) is the state with zero momentum and energy normalized to 
(xo,aj|0) = 1, and thus it is periodic too. 

In Ref. [28] this method has been applied to compute the anomalous component 
of the effective action of two- and four-dimensional fermions at finite temperature in 
the presence of arbitrary vector and axial gauge fields and scalar and pseudoscalar 
fields on the chiral circle. The computation is carried out to leading order in a 
suitable commutator which preserves chiral symmetry. As is well known, at zero 
temperature the gauge Wess-Zumino-Witten (WZW) action, which saturates the 
chiral anomaly, is the only leading contribution in the anomalous sector; further 
terms must be Lorentz and chiral invariant and they vanish identically unless they 
have more gradients thereby being sub-leading terms. At finite temperature the 
situation is different since Lorentz invariance is partially broken and this allows to 
have new chiral invariant contributions^ of the same order as the WZW action. In 
particular these terms modify the vr — > 77 amplitude which is no longer determined 
by the chiral anomaly [30]. The calculation in [28] confirms previous results [31,32] 
that this amplitude vanishes in a chiral symmetric phase (see also [33,34]). 

In Ref. [29] the same technique is applied to the study of 2 + 1-dimensional 
fermions at finite temperature in the presence of arbitrary background gauge fields. 
The use of the ^-function regularization guarantees the gauge invariance of the re- 
sult under topologically small and large transformations. This has allowed to solve 
a long standing puzzle, namely, the apparent renormalization of the Chern-Simons 
coefficient at finite temperature, which has been shown to be a perturbation theory 
artifact [35] . In Ref. [29] all ultraviolet divergent terms of the effective action, within 
a strict gradient expansion, have been computed. The result preserves gauge and 
parity symmetries (up to the standard temperature independent parity anomaly) 
and display the correct 27ri multivaluation introduced by the Chern-Simons term. 
The known exact result for massless fermions [36] is also reproduced. 

ACKNOWLEDGMENTS 

One of us (E.R.A.) acknowledges the organizers for the stimulating atmosphere 
during the workshop. This work is supported in part by funds provided by the Span- 

^) It is well established that the chiral anomaly is temperature independent (see e.g. [19]). 



ish DGICYT grant no. PB95-1204 and Junta de Andalucia grant no. FQM0225. 

REFERENCES 

1. S. Weinberg, Physica 96 A (1979) 327. 

2. See e.g. W. Marciano and H. Pagels, Phys. Rep. 3 (1978) 137 and references therein. 

3. H. Pagels and S. Stokar, Phys. Rev. D 20 (1979) 2947. 

4. J. Gasser and H. Leutwyler, Ann. Phys. NY 158 (1984) 142. 

5. For a review see e.g. C. D. Roberts and A. G. Wilhanis, Prog. Part. Nucl. Phys. 33 
(1994) and references therein. 

6. C. D. Roberts, R. T. Cahill and J. Prashifka, Ann. Phys. (N.Y.) 188 (1988) 20. 

7. B. Holdom, J. Terning and K. Verbeek, Phys. Lett. B 232 (1989) 351. 

8. R. D. Ball and G. Ripka, in Many Body Physics, World Scientific, 1994. C. Fiolhais, 
M. Fiolhais, C. Sousa and J. N. Urbano (Eds). 

9. B. Goin, W. Broniowski and G. Ripka, Phys. Lett. 437 (1998) 24. 

10. L.L. Salcedo and E. Ruiz Arriola, Ann. Phys. (N.Y.) 250 (1996) 1. 

11. E. Ruiz Arriola and L. L. Salcedo, Phys. Lett. B 450 (1999) 225. 

12. H. D. Politzer, Nucl. Phys. B 117 (1976) 397. 

13. R.T. Seeley, Amer. Math. Soc. Proc. Symp. Pure Math. 10 (1967) 288. 

14. S.W. Hawking, Comm. Math. Phys. 55 (1977) 133. 

15. E. P. Wigner, Phys. Rev. 40 (1932) 749. 

16. N. G. Pletnev and A. T. Banin, hep-ph/9811031. 

17. L. Alvarez-Gaume and P. Ginsparg, Ann. Phys. (N.Y.) 161 (1985) 423. 

18. S.L. Adler and W.A. Bardeen, Phys. Rev. 182 (1969) 1517. 

19. A. Gomez Nicola and R.F. Alvarez-Estrada, Int. J. Mod. Phys. A 9 (1994) 1423. 

20. W.A. Bardeen, Phys. Rev. 184 (1969) 1848. 

21. T.E. Clark and S.T. Love, Nucl. Phys. B 223 (1983) 135. 

22. J. Minn, J. Kim and C. Lee , Phys. Rev. D 35 (1987) 1872. 

23. J. Bijnens and J. Prades, Phys. Lett. B 320 (1993) 130. 

24. E. Ruiz Arriola and L.L. Salcedo, Nucl. Phys. A 590 (1995) 703. 

25. W.A. Bardeen and B. Zumino, Nucl. Phys. B 244 (1984) 421. 

26. H. Leutwyler, Phys. Lett. B 152 (1985) 78. 

27. E. Ruiz Arriola and L. L. Salcedo, hep-th/9811073. 

28. L.L. Salcedo, Phys. Rev. D58 (1998) 125007. 

29. L.L. Salcedo, Nucl. Phys. B549 (1999) 98. 

30. R.D. Pisarski and M. Tytgat, Phys. Rev. D56 (1997) 7077. 

31. R.D. Pisarski, Phys. Rev. Lett. 76 (1996) 3084. 

32. R. Baier, M. Dirks and O. Kober, Phys. Rev. D54 (1996) 2222. 

33. S. Gupta and S.N. Nayak, hep-ph/9702205. 

34. F. Gelis, Phys. Rev. D59 (1999) 076004. 

35. G. Dunne, K. Lee and C. Lu, Phys. Rev. Lett. 78 (1997) 3434. 

36. L. Alvarez-Gaume, S. Delia Pietra and G. Moore, Ann. Phys. 163 (1985) 288. 



